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Results  of  this  stud;  on  features  of  flow  field  and  atmosphe.ic 
diffusion  parameters  are  described  the  paperAThe  analysis  shows 
that  there  are  some  obvious  and  significant  flow  patterns  in 
mountainous  area,  such  as  valley  breeze,  channel  wind,  fumigation  and 
thermal  effect  of  industrial  sources  etc.  The  parameters  in  three 
dimensions  have  been  obtained  by  probing  on  site. 

In  this  paper,  some  of  the  problems  associated  with  the  calculation 
of  parameters  by  the  method  of  sliding  average  using  limited  number 
of  data  are  discussed.  By  using  mathematical  simulation,  the  relationship- 
between  sampling  time  and  time  interval  of  the  sliding  average  has 
been  analyzed. 
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A  STUDY  ON  FEATURES  OF  FLOW  FIELD  AND  ATMOSPHERIC  DISPERSION 
PARAMETERS  IN  MOUNTAINOUS  AREA  BY  BALANCED  BALLOONS 
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_  This  paper  introduces  the  results  of  a  study 
on  features  of  planetary  boundary  layer  flow  field 
and  rules  of  atmospheric  dispersion  at  an  industri¬ 
al  base  in  a  mountainous  area  by  using  the  balanced 
balloon  observation  method.  The  analysis  shows  that 
there  are  obvious  flow  field  features  such  as  valley 
breeze,  channel  wind,  fumigation  and  effects  of 
local  industrial  heat  sources.  Three-dimensional 
atmospheric  dispersion  parameters  have  been  obtained 
by  on-site  probing.  In  this  paper,  problems  associ¬ 
ated  with  the  calculation  of  atmospheric  dispersion 
parameters  by  the  method  of  sliding  average  using 
limited  number  of  data  are  also  discussed.  By 
using  the  method  of  mathematical  simulation,  the 
relationship  between  the  time  interval  and  sampling 
time  of  the  sliding  average  has  been  analyzed. 


I.  Preface 


The  balanced  balloon  can  be  used  as  an  indicating  particle  to 
show  the  movement  of  air  particles.  From  the  trajectory  of  its  move 
ment,  the  conditions  of  air  flow  field  can  be  understood  and  the 
atmospheric  dispersion  capability  estimated.  From  December  1981  to 
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January  1982,  we  conducted  observation  experiments  of  balanced  bal¬ 
loons  on  both  banks  of  the  Jingsha  River  at  Dukou  area  in  Sichuan 
Province  and  obtained  good  results. 

A  regular  No.  20  weather  balloon  field  with  hydrogen  and  carbon 
dioxide  at  a  ratio  of  1. 7:1.0  was  used  as  the  balanced  balloon.  The 
balloon  was  carried  by  a  captive  balloon  to  an  altitude  of  250  meters 
and  released  through  a  special  releasing  device.  The  bivane  anemom¬ 
eter  was  used  to  trace  its  movement  and  readings  were  taken  at  15 
second  intervals.  In  order  to  assure  a  balanced  mass  of  the  balanced 
balloon,  various  specialized  performance  tests  were  conducted  before 
the  experiment.  The  operation  conditions  and  observation  procedures 
were  strictly  controlled  during  the  tests,  and  the  balance  performance 
was  reliable^1^. 

II.  Track  Analysis  And  Features  Of  Mountainous  Flow  Field 

The  projection  method  and  vector  method  were  used  successively 
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to  compute  the  trajectories  of  the  balanced  balloons  .  There  was 
essentially  no  difference  in  the  balloon  location  determined  by  the 
two  methods;  the  only  difference  was  in  computation.  Actual  computa¬ 
tions  have  shown  that  the  vector  method  results  in  higher  accuracy 
and  that  it  can  directly  use  short  lines  to  indicate  index  errors. 

Its  requirements  for  setting  up  a  baseline  are  not  stringent.  The 
computation  results  can  be  markedly  improved  under  conditions  when 
the  balloon  crosses  the  baseline  or  its  trajectories  become  more 
complex.  The  vector  method  is  better  suited  for  computation  using 
a  computer,  thus  it  has  supeiority.  Certainly  its  computations  are 


more  involved  and  not  as  convenient  as  those  of  the  projection  method. 


After  obtaining  the  location  of  the  balloon,  the  coordinate  axes 
were  set  according  to  geographical  directions  with  the  releasing  point 
as  the  origin.  Then  the  horizontal ' and  vertical  movement  trajectories 
were  plotted  on  the  grid  paper  and  the  movement  distance  of  the  bal¬ 
loon  during  a  certain  period  could  be  directly  determined  from  the 
plot.  Using  the  above  method  and  a  selection  process,  135  effective 
trajectories  were  obtained  (see  Fig.  1)  for  the  use  of  the  analysis 
and  computation  of  the  atmospheric  dispersion  parameter.  Several 
typical  flow  field  features  in  a  mountainous  area  were  discovered 
after  the  analysis: 


Fig.  1.  Horizontal  trajectories  of  the  balloons  during  the  experiment 
Key:  (1)  Nong  Nong  ditch;  (2)  Nala  ditch;  '3)  Renke  Ditch; (4)  Jingsha 
River;  (5)  releasing  point. 
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1.  Features  Of  Valley  Breeze 


In  Fig.  1,  if  an  observation  point  is  used  as  the  origin  and  let 
the  x  axis  point  to  the  north  and  the-  y  axis  to  the  east,  then  the 
number  of  balloon  trajectories  in  quadrants  I,  II,  III  and  IV  of  the 
normal  coordinant  system  are  4,  49,  36  and  46,  respectively.  The 

analysis  indicates  that  the  releasing  time  of  balloons  in  quadrants 

XV 

II  and  Illyjwere  mostly  between  afternoon  and  dusk,  and  the  atmospheric 
structures  were  mostly  unstable  while  a  few  were  neutral  after  sunset. 
The  releasing  time  of  balloons  in  quadrants  IV  were  mostly  between 
morning  and  noon,  and  the  atmospheric  structures  were  mostly  in  the 
interim  period  between  stable  and  neutral.  It  can  be  observed  from 
the  daily  change  in  local  wind  patterns  that  the  westwardly  flight 
trajectories  in  quadrants  II  and  III  coincide  with  valley  breeze  and 
those  of  the  southeastwardly  in  quadrant  IV  coincide  with  the  mountain 
breeze.  Around  11  a.m.  local  time  is  exactly  the  time  when  the  night¬ 
time  mountain  breeze  starts  changing  to  valley  breeze  and  11  p.m.  is 

the  time  when  the  daytime  valley  breeze  starts  changing  to  mountain 

.  [4] 

breeze 

2.  Fumigation  Process 

Since  the  fumigation  process  often  occurs  at  the  same  time  as  the 
mountain  valley  breeze  changing  time  (around  11  a.m.)  at  the  Dukou  area, 
high  ground  level  air  pollutant  concentration  results.  This  fumigation 
process  can  be  clearly  seen  on  the  flight  trajectories  of  the  balanced 
balloons.  The  local  terrain  calls  for  the  flight  trajectories  of  this 
process  to  be  all  in  quadrant  IV  and  between  8,  9  a.m.  to  around  11  a.m. 


There  were  10  such  trajectories  during  the  experiment.  Figure  2  shows 
two  typical  groups  of  examples.  The  existence  and  development  of  a 
temperature  inversion  layer  could  be  seen  from  the  low  altitude 
temperature  probing  conducted  simultaneously.  The  trajectory  analysis 
shows:  (1)  the  occurrence  of  fumigation  process  will  certainly  bring 

air  pollutants  at  high  stratus  to  the  ground;  (2)  the  balloon  changes 
flight  direction  at  a  certain  altitude  in  the  low  stratus,  and  the 
altitude  for  direction  change  rises  as  time  goes  on;  (3)  as  time  goes 
on  the  touchdown  point  of  the  balloon  moves  closer  and  closer  to  the 
releasing  point. 


Fig.  2.  Trajectories  of  balloon  under  a  typical  fumigation  process 
No.  29:  WNW,  £=0.45  m/s,  1981.12.28,  10:55'15" 

No.  30:  WNW,  £=0.48  m/s,  1981.12.28,  1 1 : 27 • 1 5" 


3.  Effects  Of  Local  Industrial  Heat  Sources 


The  two  groups  of  flight  trajectories  in  Fig.  3  show  that,  after 
releasing,  the  balloon  started  to  descend  due  to  the  effects  of  down¬ 
ward  air  flow.  But  when  it  flew  over  the  high  temperature  plume  of  a 
large  industrial  source,  it  was  being  pushed  upward  and  within  1  to  2 
minutes,  it  rose  155  meters  and  70  meters,  respectively.  The  atmo¬ 
spheric  structure  was  stable  at  this  time.  Thus,  we  used  the  plume 
rise  formula  under  corresponding  conditions  to  calculate  plume  rises, 


SW’*? 


V$! 

|« 

m 


and  the  values  were  123  meters  and  78  meters,  respectively.  These 
were  very  consistent  with  the  heights  of  balloon  rise  at  that  location. 
This  phenomenon  indicates,  at  the  very  least,  that  the  effects  of  local 
large  industrial  heat  sources  cannot  be  ignored. 


Fig.  3.  Effects  of  industrial  heat  sources 


4.  Effects  Of  Channel  Wind 


It  can  be  seen  from  the  two  examples  shown  in  Fig.  4  that  after 
the  balloons  crossed  the  Jingsha  River  from  east  to  west,  their  flight 
directions  were  consistent  with  the  channel  direction  except  that  the 
latter  half  of  the  flight  of  No.  77  deviated  slightly.  This  was  be¬ 
cause  the  balloon  kept  rising  as  it  flew  over  the  river.  When  it 
reached  above  450  meters  it  started  to  deviate  from  the  channel  direc¬ 
tion.  This  shows  that  the  effects  of  the  channel  can  reach  an  altitude 
between  400  and  500  meters.  The  effects  of  channel  downward  air  flow 
were  slightly  visible.  But  since  the  atmospheric  conditions  during 
the  flight  of  the  two  groups  of  trajectory  were  quite  unstable,  the 
effects  of  air  flow  were  greater  and  part  of  the  effects  of  downward 
air  flow  were  covered  up. 


Fig.  4.  Effects  of  Jingsha  River  valley 
Key:  (1)  Jingsha  River. 

III.  Determination  Of  Atmospheric  Dispersion  Parameters 

1.  Discussion  Of  Computation  Methods 

The  Lekkau  wind  velocity  pulse  data  were  obtained  from  the  flight 
trajectory  of  a  single  balloon.  Based  on  the  Taylor  theorem,  the 
three-dimensional  atmospheric  dispersion  parameters  are  computed  from 
the  following  formulae: 

°l=  («'*)*.  ,Tl 

a}  =(«'•),  jT'  (2) 

where  u',  v'  and  w1  are  the  wind  velocity  pulses,  respectively.  The 

subscript  T  represents  the  sampling  time,  T  is  the  release  (or  travel) 
-  2  -  ?  _  2 

time.  Here  (u'  )  T,  (v'“)^  T  and  ( w '  )  ^  T  represent  the  square 

deviations  of  sliding  average  velocity ,  respectively ;  that  is  to  cal¬ 
culate  the  sliding  average  of  wind  velocity  data  over  T  with  respect 
to  T,  then  the  square  deviations  are  obtained,  respectively. 

When  using  a  limited  number  of  wind  velocity  pulse  sliding  aver¬ 
ages  to  obtain  the  dispersion  parameters,  their  calculated  values  have 


a  tendency  to  drop  irregularly  as  the  time  interval  reaches  a  certain 


value.  As  the  interval  continues  to  increase,  a  decrease  even  in  the 
absolute  values  of  the  calculated  values  will  result.  Therefore,  it 
has  been  pointed. out  by  people^ ^  that  the  maximum  travel  time  should 
only  be  selected  as  10-20%  of  the  sampling  time.  We  believe  that, 
from  the  energy  spectrum  angle,  this  suggestion  is  reasonable.  There 
are  cutting  effects  on  the  low  and  high  frequencies  of  the  energy 
spectrum  by  the  sampling  time  with  a  limited  length  and  average  time 
with  a  fixed  interval.  We  computed  the  effects  of  limited  sampling 
time  (TV  oo)  and  average  time  (T^O)  on  square  deviation  and  found 

9  2 

that:  only  when  TV  5  can  the  closeness  of  rj  m  and  d  “  _ 

w  T  •  T  u  oo  .  0 

be  assured.  Therefore,  it  is  appropriate  to  select  the  maximum  slid¬ 
ing  average  interval  T  as  10-20%  of  'V . 

In  order  to  find  the  relationship  between  the  sliding  time  inter¬ 
val  and  sampling  time  corresponding  to  the  above  phenomenon,  the 
mathematical  simulation  method  was  used  to  conduct  further  analysis. 
Assume  there  is  a  vector  V  with  unit  length  and  it  changes  direction 
with  respect  to  t.  It  changes  randomly  within  a  360°  range  in  space. 
If  the  y  direction  component  of  this  vector  at  any  time  in  an  unit 
time  interval  is  selected  as  the  simulated  y  direction  wind  velocity 
pulse,  then  its  standard  deviation  can  be  obtained  by  using  the  slid¬ 
ing  average  method.  If  the  vector  moves  an  unit  distance  in  the  x 
direction  within  an  unit  time,  and  when  the  samples  are  selected  large 
enough,  then  it  can  be  proved  that: 

<3, 
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Twenty  groups  of  the  above  number  were  generated  randomly  by  computer 
with  each  group  consisting  of  50  and  80  random  numbers.  The  standard 
deviation  d  of  different  sliding  time  was  obtained  by  using  the 
sliding  average.  The  calculated  resu-lts  show  the  following  conclusion 
(1)  the  calculated  standard  deviations  of  either  50  or  80  numbers  are 
smaller  than  those  calculated  from  formula  (3).  The  larger  the  slidin 
time  interval,  the  greater  the  decrease;  (2)  the  relative  velocities, 
which  decrease  as  the  time  interval  increases,  of  two  numbers  are 
in  fair  agreement.  As  they  reach  10%,  20%  and  30%  of  each  individual 
length,  the  average  values  decrease  by  4.0%,  5.1%;  14.2?6,  13.7%  and 
24.5%,  25.0%  respectively;  (3)  when  sliding  to  30%  (for  50  numbers) 
and  40%  (for  80  numbers)  of  the  total  length,  the  absolute  values  of 
the  average  values  start  to  drop;  (4)  the  standard  deviation  for  each 
group  is  obtained  successively  and  it  is  found  that  there  is  greater 
difference  between  each  group.  When  reaching  20%  of  the  total  length, 
the  maximum  values  are  about  3  times  larger  than  the  minimum  values. 
Figure  5  shows  these  effects. 


Fig.  5.  Effect  limiting  number  on  sliding  average 

To  use  what  method  in  computing  the  average  values  of  d  under 
different  atmospheric  conditions  is  also  an  important  question.  We 
tried  six  average  methods  and  compared  the  results,  and  the  method 
with  the  best  results  was  selected;  i.e.. 


first  regrouped  and  converted 


the  data  to  obtain  the  coefficient  A  and  index  b  in  the  index  form  of 
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dispersion  parameter,  then  computed  the  standard  deviation  of  each 
group  according  to  its  index  form,  and  lastly  took  the  average  value 
and  converted  it  to  obtain  new  A  and  b  values. 


2.  Results 


The  atmospheric  dispersion  parameters  were  computed  under  five 
categories:  stable  (W)  ,  nuetral  (W) ,  neutral  (E),  unstable  (E)  and 

breeze  which  were  divided  according  to  atmospheric  conditions,  wind 
direction  and  wind  velocity.  The  categorization  of  atmospheric  condi¬ 
tions  is  mainly  based  on  data  from  on-site  low  altitude  temperature 
probings  and  determined  by  the  following  standards:  i.e., 
C^f]Z<-°s'c/To0"  is  unstable;  -0.5t/lD0m<f"  <0.4t/100m  is  stable 


Neutral  (W)  includes  a  weak  stable  category  and  neutral  (E)  includes 
a  weak  unstable  category.  The  calculated  results  of  dispersion 
parameters  for  each  category  are  shown  in  Table  1.  It  can  be  seen 


from  Table  1  that:  (1)  under  breeze,  <S  x  and  &  are  closer  (value 


and  trend  of  change),  and  they  are  both  3  times  larger  than  Cf ; 


(2)  under  breeze,  both  O'  and  ( T  are  much  larger  (3-5  times)  than 

x  y 


those  under  unstable,  neutral  and  stable  categories.  (T  is  within 


the  same  magnitude  as  that  under  unstable  category;  (3)  (j  x,  (5  and 


(S’  under  the  unstable  category  are  larger  than  those  under  the  neutral-! 


and  stable  categories.  For  the  unstable  category  itself,  Cf  and  (j 

y  z 


are  a  little  larger  than  (T  ,  and  (S'  and  CT  are  about  the  same 

x  y  z 


magnitude;  (4)  6"  and  (S'  of  the  stable  category  are  larger  than  those 

y  z 


under  the  neutral  category,  6  is  about  the  same  magnitude  as  those 
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under  the  neutral  category;  (5)  for  the  neutral  category  itself, 


(f  v  >  <fx  >  (f  z*  For  the  stable  category,  Cfy  >  <TX  and  z>  and  Cf 


and  <5  are  about  the  same  magnitude.  These  trends  are  in  fair 


agreement  with  the  actual  conditions  in  the  Dukou  area.  The  breeze  at 


Dukou  usually  occurs  during  the  changing  period  of  the  mountain  and 


valley  breeze.  At  such  time,  the  wind  direction  changes  erratically 


and  the  balloon  trajectories  appear  to  curve  horizontally.  This  is 


the  primary  reason  that  results  in  large  breeze  dispersion  parameters. 


Except  that  the  dispersion  parameters  under  unstable  category  show 


a  trend  of  equal  characters  in  all  directions,  this  trend  is  not 


obvious  for  those  under  other  categories.  This  situation  is  reasonable 


for  the  Dukou  area  where  the  mountain  is  high,  the  valley  is  deep  and 


the  wind  velocity  at  high  altitude  is  small.  CT ^  and  CT ^  of  the 


stable  category  are  larger  than  those  under  the  neutral  category 


unexpectedly.  It  was  discovered  after  inspecting  the  data  of  each 


group  under  stable  category  that;  the  releasing  altitudes  of  the 


balloons  in  these  groups  were  mostly  located  in  the  successively 


declining  layer  or  the  isothermal  layer  of  the  mountainous  multi-layer 


temperature  inversion  structure,  i.e.,  they  were  still  flying  in  the 


unstable  or  neutral  layer  where  there  was  string  temperature  inversion 


layer  in  the  upper  portion  and  either  no  temperature  inversion  layer 


or  a  thin  temperature  inversion  layer  in  the  lower  portion.  Thus, 


the  temperature  inversion  layer  only  limited  the  balloon's  ascent, 


yet  did  not  exert  very  much  influence  on  the  balloon's  descent. 


Since  the  influence  was  not  uniform,  the  balloon  was  actually  flying 


most  of  the  time  in  the  neutral  or  unstable  category  causing  the 


dispersion  parameters  to  be  still  quite  large. 
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Uable  1.  Dispersion  parameters  under  various  conditions  f  Q  =  Ax 
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Fig.  6.  Comparison  of  curves 
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— x —  Wind  tunnel  simulation  at  Dukou  — A —  Photographed  at  Dukou 
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the  absolute  values  of  dispersion  parameter  start  to  drop.  This 


requires  that  ue  obtain  as  long  a  sampling  time  as  possible.  Under 
conditions  where  the  actual  sampling  time  is  long  enough,  it  is  better 
to  take  the  sliding  average  time  interval  to  be  10%  of  the  total  sam¬ 
pling  time,  but  not  to  exceed  20%  at  the  most;  otherwise  there  will 
be  larger  errors  between  the  calculated  dispersion  parameter  values 
and  the  observed  values.  When  conducting  balanced  balloon  observa¬ 
tions  in  a  mountainous  area,  the  observation  might  be  terminated 
prematurely  due  to  the  obstruction  by  terrain  or  land  objects  causing 
the  data  collected  to  be  insufficient.  Therefore,  the  rate  of 
selective  elimination  of  the  data  group  tends  to  be  higher.  This  is 
a  problem  which  demands  serious  attention  and  improvement. 

The  study  in  this  paper  was  part  of  the  Dukou  environmental 
quality  overall  evaluation  research  subject.  During  the  course  of 
this  study,  full  support  and  assistance  were  received  from  the 
coordination  group.  The  Bureau  of  Meteorology  of  Dukou  City  not  only 
provided  us  with  part  of  the  high  quality  meteorological  data,  but 
also  supported  and  participated  in  part  of  the  field  work.  We  hereby 
express  our  special  thanks  to  all  of  them. 
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